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Abstract: A highly excited Friedmann universe filled with a scalar field and radia- 
tion has been considered. On the basis of a direct solution to the quantum-mechanical 
problem with a well-defined time variable, it has been shown that such a universe can 
have features (energy density, scale factor, Hubble constant, density parameter, matter 
mass, equivalent number of baryons, age, dimensions of large-scale fluctuations, ampli- 
tude of fluctuations of cosmic microwave background radiation temperature) identical to 
D ■ those of the currently observed Universe. 
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1. Available cosmological data suggest that, from the point of view of quantum the- 
ory, the currently observed Universe is likely to be in a highly excited state [1]. This 
^ ■ is confirmed by estimates of the number of the quantum state that corresponds to its 

averaged motion as a discrete unit [1-4]. In view of this, it is important to investigate cos- 
mological systems featuring gravitational and matter fields and occurring in states with 
large quantum numbers by proceeding from a direct solution to the relevant quantum- 
mechanical problem. 

A model of the Friedmann universe filled with a uniform scalar field has been proposed 
in [4]. This model, which is appropriate for constructing a quantum theory, features a 
well-defined time variable. The reference frame was specified there with the aid of a 
subsidiary matter source in the form of radiation (relativistic matter of any nature) that 
was assumed to be initially present in the cosmological system along with a scalar field 
that forms a nonzero cosmological constant in the early universe. The evolution of the 
universe filled with not only scalar field but also with radiation differs from that which is 
realized in the absence of radiation. The main difference lies here in the emergence of a 
new region that is accessible to a classical motion and which is bounded by the potential 
barrier existing in the system of scalar and gravitational fields. A quantum universe 
involving a slowly varying scalar field and occurring in low-lying (quasistationary) states 
has been analyzed in [4], where it has been shown that the dynamical model proposed 
there is compatible with the currently prevailing ideas of the early Universe. In this 
study, we will consider a quantum Friedmann universe with large quantum numbers 
characterizing the possible physical states of the gravitational and matter fields involved. 
On the basis of a solution to the quantum mechanical problem, it is shown that the 
universe in highly excited states can have features (energy density, scale factor, Hubble 
constant, density parameter, matter mass, equivalent number of baryons) identical to 
those in the currently observed Universe. The problem of the age of the Universe and 
the possible new mechanism that could generate fluctuations of the metric due to a finite 
width of quasistationary states and to the anisotropy of cosmic microwave background 
radiation are discussed in the Appendices. 
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2. The wave function of the quantum Friedmann universe filled with radiation and a 
uniform scalar field <fi specified by the potential V(<j>) is determined by the Schrodinger 
equation [4] 

2ifr*=(%-l8i-U\*, (1) 

where T is a privileged time coordinate related to the synchronous proper time t by the 
coordinate condition dT/dt = 1/a, a being the scale factor, while 

U = a 2 -a 4 V{<j>) (2) 

plays the role of the effective interaction potential in the cosmological system being 
considered. Here, we use the system of units in which I = ^2Gh/3Tcc 3 = 1 and (f) = 
^/3c 4 /87rG = 1.Q 

Possible solutions to equation (1) are determined by the properties of the scalar field 
involved. If V(<p) is everywhere positive definite, we can see that, in the variable a, 
the potential U has the form of a barrier with height U max = and width Aa = 
(12(E) — a\(E), where a±(E) and (12(E), a%(E) < (12(E), are the classical turning points 
determined from the condition U = E. With the aid of the equation of motion for the field 
0, it can be shown that, if the field satisfies the condition |<9 f 2 0| -C \^\, the contribution 

of the operator -^d 2 can be approximated by the term — ^ > wnere H = ^ is 

the Hubble constant; that is, this contribution is equivalent to the squared addition to the 
interaction U of the gradient of the potential V(<j)). In this effective potential, stationary 
states cannot exist in the region a < a\. If, however, V(<fi) 1, quasistationary states 
with lifetimes exceeding the Planck time can exist within the barrier [4]. The positions 
and widths of such states are determined by the solutions to equation (1) that satisfy 
the boundary condition in the form of a wave traveling toward greater values of a. 
A general solution to equation (1) can be represented in the integral form 



*(a,0,r)= / dEe^ 1 C(E)^ E (a,4>), (3) 
Jo 

where the function C(E) characterizes the E distribution of the states of the universe 
at the instant T = 0, while ^(a,^) and E are, respectively, the eigenf unctions and the 
eigenvalues for the equation 

d 2 a + ^dl + U-E\^ E = 0. (4) 



3. Let us now consider a quantum universe where 
(4) can then be represented as 



1 dV 

v d<p 



2 

C 1. A solution to equation 



ip E (a,^)= deip e (a,(j)) f e (<f>;E), (5) 







1 In order to go over to the system of units where K = c = 1, we must use the relation lip 2 = 
y/3/32TT 3 m p for the energy and the relation (4>/l) 2 — (9/16) for the energy density, m p being the 
Planck mass. 
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where (p e and e are, respectively, the eigenfunctions and eigenvalues for the operator 
[— d 2 + U] of the adiabatic approximation that correspond to continuum states at a fixed 
value of the field 0. The functions tp e can be normalized to the delta function 8{e — e'). 
Their form greatly depends on the value of e. For e < U max , there are quasistationary 
states with e = e n = e n + i T n in the system [4], and the main contribution to the integral 
in (5) over the interval < e < U max comes from the values e ~ e n and a < R, where 
R > a 2 (e). Here, the wave function has the form 

<p e = A(e)<p?\ (6) 

where the function A(e) has a pole in the complex plane of e at e = e n , while is 
the solution over the interval < a < R that is regular at the point a = 0, normalized 
to unity, and weakly dependent on e. Proceeding in a way similar to that adopted in 
the theory of quasistationary states for short-range potentials [5], we can show that the 
quantity |A(e)| 2 can be approximated by the delta function <5(e — e n ) in the case of the 
potential (2) as well. In this approximation, expression (5) can eventually be reduced to 
the expansion 

POO 

ipE(a, 0) = ^2 fn(a, 0) / n (0; E) + / de(p € (a,(p) f e (<f>;E), (7) 

^ J U max 

where <p n = (p^J for < a < R and tp n — for a > R (a state of this type was considered 
in [6]), while f n {<t>] E) = da f^(a, 0) iPe(o>, <P)- In the limit of an impenetrable barrier, 
the function ip n reduces to the wave function of a stationary state with a definite value 
of e n . In the case of V 1 considered here, the contribution of the integral to the 
expansion in (7) can be disregarded, and the quantities f n ((p;E) can be interpreted as 
the amplitudes of the probability that the universe is in the state f n {4>', E) with a given 
value of the field 0. They satisfy the set of differential equations 

dlf n + J2 K nn>(<f>;E)f n ,=0, (8) 
n' 

where 

K nn ' = ((p n \dl\(p n >) + 2 (<p„|fy|<p„')fy + ^ ((p n \d 2 | w) (^n' - E) . (9) 

Over the time interval AT < we can disregard the possibility of decay and consider the 
quasistationary state as a stationary state that arises in place of the quasistationary state 
when the decay probability tends to zero. The wave function ip n of such a stationary state 
and the corresponding eigenvalue e n can be found by perturbation theory by considering 
the interaction a 4 V(0) as a small perturbation against a 2 (in the region a < ai, we have 
a 2 V < 1). This yields 

= \n)~j \ VN(N -1){N- 2){N - 3) \n - 2) + 
+ y/N(N-l) (N-^j |n-l>- 

- ^/(N +1)(N + 2) (n + ^ \n + 1) - 

- I V(N + 1)(N + 2)(N + 3)(N + 4) \n + 2) + O (V 2 ) , (10) 
8 
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= el - \v [2N(N + 1) + 1] 
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( +0{V 3 ), (11) 

where A" = 2n + 1, while e° = 2iV + 1 and |n) are, respectively, an eigenvalue and 
the corresponding eigenfunction of the operator [— d 2 + a 2 ] . Quasistationary states are 
realized for V < 0.08 = 4.5 x 10~ 2 mp and are characterized by values e n > 2.6 and 
r n <C 0.3 [4]. The universe can undergo a tunnel transition to the region a > a<i (e n ) from 
any quasistationary state. 

In the early universe, the quantity V((f>) specifies the vacuum energy density, which 
determines the cosmological constant A at that era [7, 8]. In our Universe, the cosmo- 
logical constant is very small. In the model featuring matter in the form of one scalar 
field, the reduction of the cosmological term can be described in terms of the potential 
V(4>(t)), which decreases with time [3, 9]. A similar behavior of V((f>(t)) is suggested by 
the results of investigations within inflationary models [8, 9]. It was shown in [4] that, 
with a nonzero probability, a quantum universe initially filled with radiation and a scalar 
field whose potential V(<f>(t)) decreases with time can evolve in the region bounded by 
the barrier. The expansion occurs here owing to transitions from lower to higher states 
via the interactions of the scalar and gravitational fields involved. In the zeroth ap- 
proximation, the evolution of the universe can be described by considering transitions 
between unperturbed states \n) and \n + 1) under the effect of the interaction a 4 V [4]. 
In a more rigorous approach that takes into account variations in V(<f>), the universe 
expands through transitions between the states f n and (p n +i due to the gradient of the 
potential V(<f>). 

As the potential V(<f>) decreases, the number of quantum states in which the universe 
can occur increases. As a result, the competition between the tunneling processes and 
transitions from one state to another arises in the universe that has not had time to 
undergo tunneling. Since the decay probability decreases exponentially with decreasing 
V, the probability that the universe is excited to states with large quantum numbers n 
in the course of time is nonzero. The probability that, in course of time, the universe in 
large-n states will occur in the region outside the barrier is negligibly small, because small 
values of V correspond to such states. In the limit V — >■ 0, the universe is completely 
locked in the region within the barrier with e n tending e°. 

4. Let us investigate the properties of the universe in states with n 1. In such 
states, |V| 1 and T n ~ 0, and we can use the wave function (10) in calculating the 
quantities K nn i in (9). As before, we assume that the potential V(4>) changes slowly in 

these states as well ( ^ « 1). The derivatives d^ipn in (9) can then be disregarded. 
Taking this into account and going over to the limit n 1, we can then reduce equation 
(8) to the simplified form 

&f n + ul(<f>;E)f n = 0, (12) 



where 

LU 2 n (<P;E) = 2N 2 



l-^-2W ( 0) 



(13) 



The potential of the field will be chosen in the form V((f>) = ^0 2 . From the condition 
|V| 1, it follows that the mass of the fields must be constrained by the condition 
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m ^ T^f ' ^ n ^is case, only the region |0| m p corresponds to extremely large masses 
m > m p . For masses satisfying the inequality m <C m p , the field <fi can be determined over 
the entire interval < \<p\ < oo; in this case, equation (12) reduces to the Schrddinger 
equation for a harmonic oscillator. Under such conditions, the field <fi can oscillate near 
the minimum of the potential V((f>), causing the production of particles.^ Substituting the 
explicit expression for V((j)) into (13) and introducing the new variable x = (2m 2 iV 3 ) 1//4 0, 
we arrive at the equation 

d 2 J n + (z- x 2 )f n = 0, (14) 

where z = (l — ^). Equation (14) has solutions at z = 2s + l, where s — 0, 1, 2, 

Suppose that the universe occurs in a state with a definite quantum number s. The 
function /„ can then be represented as 

1 / 1 \ 1/4 

f ns = -= {Bl) s f n0 , B n f n0 = 0, f n0 = ( - J e"* 2 / 2 , (15) 

where f n o is the state vector of the universe in the nth state at s — 0, while = 
-^(x — d x ) and B n = -j^(x + d x ) are the operators that, respectively, create and anni- 
hilate particles in this state and which satisfy the conventional commutation relations 
[B n , £>t] = 1 and [B n , B n ] = [B^, B^} = 0. Since n ^> 1, small changes in n do not affect 
the physical state of the universe where there are s particles. 
The condition of quantization of E has the form 

E = 2N - (2N) l/2 (2s + l)m. (16) 

For m m p and small s, the presence of the field has no effect on the properties 
of the universe because, in this case, the approximate equality E ~ 2N holds to a high 
precision, so that the universe is dominated by radiation. A transition from the radiation- 
dominated universe to the universe where matter (in the form of particles produced by 
the field 0) prevails occurs when, owing to an increase in the number of particles, the 
second term in (16) becomes commensurate with the first one. Let us consider the case 
where there are many particles in a highly excited nth state of the universe. For n ^> 1 
and s>l, the wave function has the form 

i> ns (a, 4>) = <p n {a)fns(<f>), (17) 

where 

/0\l/4 / AT^\ 

(18) 




fns(4>) = 




2 At E = 0, a similar mechanism leads to the production of particles by the inflaton field, which 
is identified with the scalar field <f> [8]- In the case considered here, the field for n ^> 1 states can be 
treated as an effective field obtained upon averaging over the internal degrees of freedom of real physical 
fields. 
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This wave function corresponds to the semiclassical solution to equation (4) and is 
normalized to unity with allowance for the fact that the probability of finding the universe 
in the region a > a<i is negligibly small. 

Let us now calculate the energy density p to t for matter and radiation in the universe 
described by the wave function (17). The energy density for a classical field and 
radiation is determined by the Einstein equation for the ([]) component and is given by 
[4,10] 



2 , „ E 

a: 



A + V + ~v (2°) 



where is the momentum canonically conjugate to the variable 0. Assuming that all 
the quantities in (20) are operator- valued, we set ptot = (p), where averaging is performed 
with the wave function (17). Disregarding the variances (a 2 ) — (a) 2 and (a 6 ) — (a) 6 , we 
obtain 

P« = ^«> + W + ^ (21) 

Calculating the expectation values in (21), we arrive at the total energy density in the 
universe in the form of the sum of the energy densities for matter and radiation as in the 
general theory of relativity [10]: 

193 M E 

ptot ~l2 W + W { ] 

Here, = ms, and (a) = is the scale factor in the universe occurring in an n > 1 
state. The constant E and the total mass M<p are related by the equation 

E = 4(a) [(a) - M ] . (23) 

If 

(a) = M , (24) 
there is no radiation; we then have p to t = Psub and 

/ x A93 i y /2 



5. Suppose that the system being considered occurs in an n > 1, s ^> 1 state and 
that it is characterized by the scale-factor value of (a) ~ 10 28 cm. The condition in (24) 
will then be satisfied at ~ 10 56 g. This value coincides with the mass of matter in the 
observed part of our Universe [10]. On the other hand, we set to the density of cosmic 
microwave background radiation energy in the present era, p° = 2.6 x 10~ 10 GeV/cm 3 
[11]. From (23), we then find that, for (a) = 10 28 cm, we have § = l- 0.7x 10~ 5 ; 
that is, the equality in (24) holds to a high precision in this case as well. We note that, 
although the absolute value of E is not small in the present era (E ~ 10 117 [4]), it is small 
in relation to (a) 2 ~ 10 122 (all estimates are presented here in the system of units where 
/ = = 1) and can be disregarded. It is the ratio of these quantities that determines the 
accuracy to which relation (24) holds. 
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The current values of the scale factor ((a) ~ ao ~ 10 28 cm) and of the mean energy 
density (p su b ~ po ~ 10~ 5 GeV/cm 3 ) satisfy relation (25). Associating the well-known 

equality a = (n^rr^J with relation (25), we find that the effective value of the 
density parameter is Qq = 1.066; that is, the geometry of the universe with the above 
features is close to Euclidean geometry. It is characterized by the quantum-number values 
of n ~ (a) 2 ~ 10 122 and s ~ Taking the proton mass for m, we obtain s ~ 10 80 . The 
above value of n complies with existing estimates for our Universe (see [1-4]), while s is 
equal to the equivalent number of baryons [10]. The value found for Qq corresponds to 
the Hubble constant value of }{^ heory ~ 94 km/(s Mpc). This value of Hubble constant 
lies within the limits of experimental uncertainty of the value of Hq XP = 80 ± 17 km/(s 
Mpc), which was obtained with the aid of the Hubble cosmic telescope [12]. 

The approach developed here makes it possible to obtain realistic estimates for the 
age of the Universe (see Appendix 1), for the proper dimensions of the nonhomogeneities 
of matter which are consistent with those of the observed large-scale structure of the 
Universe, and for the amplitude of the fluctuations of the cosmic microwave background 
radiation temperature. The resulting estimate for the last quantity is close to the value 
extracted from experimental data (see Appendix 2). 

6. The above numerical estimates of the parameters of the universe are of an illustra- 
tive character. If we nevertheless associate them with our Universe, it can be concluded 
that the values observed in the present era for the scale factor, the mass of matter in the 
Universe, and the density of the cosmic microwave background radiation energy satisfy 
relation (23). The zero-order approximation, which corresponds to setting 7r 2 = in 
(20), leads to a very small value of fio — 0.1, but it is very close to the lower boundary of 
the uncertainty interval for Qq [H]- Although the potential V(<p) undergoes only small 
variations in response to changes in the field 0, the field <fi itself changes fast, oscillating 
about the point = 0, so that the approximation in which 7r 2 = is invalid. The applica- 
tion of the present model in this approximation would result in the radiation-dominated 
universe; that is, it would not feature a mechanism capable of filling it with matter upon 
a slow descent of the potential V(4>) to the equilibrium position, which corresponds to 
the true vacuum. In states of the universe that are characterized by large values of the 
quantum numbers, the kinetic term of the scalar field ensures the density parameter value 
close to unity. 

Replacement of the entire set of actually existing massive fields by some averaged 
massive scalar field seems physically justified for states of the universe that have large 
values of n. We can see that, by and large, such an averaged field describes correctly the 
global features of our Universe. It effectively includes visible baryon matter and dark 
matter, which is globally manifested on cosmological scales via gravitational interaction. 
The status of the field <f> changes as we go over from one stage of universe evolution to 
another. In the early universe, the field ensures a nonzero value of the vacuum-energy 
density (cosmological constant) due to V(<p) values at which the equation for (p e (a,<f)) 
admits nontrivial solutions in the form of quasistationary states. In a later era, when the 
field <fi descends to a minimum of the potential V(4>) and begins to oscillate about this 
minimum, it appears to be a source of some averaged matter filling the visible volume of 
the universe, which has linear dimensions on the order of ~ (a). 

APPENDIX 1 

The time 4- corresponding to the theoretical Hubble constant value uH ieory [ s equal 
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to 10 10 yr, and the age of the Universe calculated by the standard formula t = is 
t w 7 x 10 9 yr. This result is smaller than the expected value of to — (10 -j- 20) x 10 9 
yr [11], but it is close to that which corresponds to H 3xp , t ~ 8 x 10 9 yr, highlighting 
the problem of the age of the Universe. It should be borne in mind, however, that the 
age of the universe is calculated by the formula t — jj, where q — \ (for the equation of 
state p — | or | (for p — 0), which implies that a = bt q , where the proportionality factor 
b is independent of t. In the model of the universe filled with matter [as represented by 
the field 4>(t)] and radiation, the factor b depends on the variable 4>(t), which changes 
with time. The inclusion of this dependence leads to a greater value of t - By way of 
example, we indicate that, for 2\/V t <C 1, the scale factor varies with t according to 
the law a(t) ~ (2 v /et) 1/ ' 2 , where e = e(0(t)) [4]. From the above, it follows that, in 
order to calculate the age of the universe, we must consider the transcendental equation 
t = [2H — dt lne] -1 . Since e > and since it increases with time, we have e > 0. 

In order to perform a numerical estimation, we assume that, in the time interval 
being considered, e grows in proportion to a power of time: e ~ t a . The age of the 
universe is then given by the expression t = ^jjp Since we have e ~ 10 117 and t ~ 10 61 
(in l/c units), the exponent a is a ~ 1.9 [a similar power-law dependence follows from 
the relation a 2 ~ n ~ e and from the classical expression for a(t)], whence we find that, 
if the Hubble constant set to its theoretical value HQ heory , the age of the universe is 
t — 15.2 x 10 9 yr. It should be noted that the above power-law dependence of e on t 
ensures a correct relationship between the current value of the scale factor ao and the 
age t of the universe: a ~ to ~ 10 61 - 

APPENDIX 2 



A quasistationary state with a small, but finite value of the width T does not possess a 
definite value of e. This uncertainty we denote it by Se can serve as source of fluctuations 
of the metric. Let us demonstrate this explicitly. By associating e + Se with the scale 
factor a + Sa and by using the solution to the Einstein equation in the region a < a\ from 
[4], we find that the amplitude of fluctuations of the scale factor can be represented as 

Sa 1 Se/e , . „ s 

— = 1= 1= ■ (-4-1) 

a 4 1 - tanh v / Ut/2 v / ^ 

Since Se < T, the fluctuations Sa that were generated at the early stage of the evolution 
of the Universe will take the greatest values. In order to estimate them, we adopt the 
values of V — 0.08, e = 2.6, and Se < 0.3, corresponding to the time t ~ 1 [4]. We 
then have ^ < 0.04. Since the dimension of large-scale fluctuations changed in direct 
proportion a(t), this relation has remained valid up to the present time [13]. Taking this 
into account, we find that Sa < 130 Mpc for the current value of a ~ 10 28 cm. On the 
order of magnitude, the above value corresponds to the scale of superclusters of galaxies 
[9]. Smaller values of Se are peculiar to quantum states with smaller V. The fluctuations 
Sa corresponding to them are smaller than those presented above and are expected to 
manifest themselves against the background of the large-scale structure. They can be 
associated with clusters of galaxies, galaxies themselves, and clusters of stars. Thus, the 
conclusions drawn on the basis of the model considered here are in line with the generally 
accepted concept that galaxies, their clusters, and other structures in the Universe are 
macroscopic manifestations of quantum fluctuations that have grown considerably [9]. 
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Let us estimate the amplitude of fluctuations of the cosmic microwave background 
radiation temperature, ST/T. By using the relation e = p 7 a 4 , where p 7 = (7r 2 /15)T 4 is 
the density of the cosmic microwave background radiation energy, we obtain 

ST 1 Se 5a , . . 

"17" ' < A2 > 

For VVt < 1, it follows from (A.l) and (A.2) that 

ST t 5a 

T a ivE7' (A3) 

For the time t ~ 10 5 yr ( t ~ 10 56 in //c units), which corresponds to the recombination 
of primary plasma (separation of radiation from matter), it can be found that, for the 
observed value of e = 5.2 x 10 117 (in which case W t < ^ ~ 0.7 x 10~ 3 ), the sought 

amplitude of fluctuations of cosmic microwave background temperature radiation at ^ < 
0.04 can be estimated as 

y<2.8x 10- 5 . (A4) 

Upon recombination, the fluctuations of the temperature undergo no changes; therefore, 
measurement of the quantity ST/T for the present era furnishes information about the 
Universe at the instant of last interaction of radiation with matter. The estimate in (A. 4) 
is in good agreement with experimental data on cosmic microwave background radiation 
(see [14]). 
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